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ABSTRACT 


The  problem  of  the  far  field  acoustic  wave  produced  in 
N lossless  plane  parallel  liquid  layers  by  point  source 
excitation  is  solved  in  a form  permitting  of  practical  compu- 
tations Formal  solutions  of  essentially  the  same  problem, 
but  not  adaptable  to  convenient  calculation,  have  been  ob- 
tained at  least  as  early  as  194-7  (Reference  (b))« 

Pekeris  has  solved  this  problem  in  a computable  form  for 
the  special  cases  of  two  and.  three  layers  (Reference  (a)  ), 
terming  the  type  of  solution  obtained  "the  normal  mods  fora". 
The  solution  to  be  given  here  for  the  N layer  case  is  simi- 
larly in  the  normal  mode  form,  reducing  to  the  results  obtained 
by  Pekeris  for  the  two  and  three  layer  cases. 

While  the  present  solution  is  a generalization  of  the 
form  obtained  by  Pekeris,  the  method  of  derivation  is  basically 
different.  In  particular,  a fundamentally  different  method  of 
obtaining  the  mode  excitations  is  employed.  Use  of  this 
method  not  only  greatly  simplifies  finding  the  Green's  function 
in  the  present  problem,  but  it  appears  capable  of  facilitating 
the  solution  of  various  other  excitation  problems. 


INTRODUCTION 


The  problem  of  the  acoustic  wave  produced  in  N lossless 
plane  parallel  liquid  layers  by  point  source  excitation  has 
(in  essence)  been  solved  formally  by  Lurye  in  194 7 (reference 
(b)  ) . However  the  first  solutions  of  this  problem  which  were 
in  a form  permitting  of  practical  computation  were  obtained 
by  Pekeris  in  1948  for  the  special  cases  of  two  and  three 
layers  (reference  (a)).  Ytfe  will  obtain  here,  by  basically 
different  techniques,  a generalization  of  Pekeris'  form  of 
solution  to  the  N layer  case. 

Each  of  the  parallel  layers  is,  in  itself,  uniform  and 
characterized  by  a density  /?  and  wave  velocity  Cy  . What 
will  be  called  the  top  layer  is  bounded  by  an  acoustic  free 
surface  while  the  bottom  layer  is  a liquid  of  infinite  depth 
(forms  a half  space).  We  will  consider  the  particular  case 
where  the  source  is  in  the  top  layer,  the  technique  being 
essentially  the  same  for  the  source  in  any  layer.  (In  fact, 
if  the  solution  is  known  for  the  source  in  the  top  layer,  it 
can  immediately  be  found  for  the  source  in  any  layer  with  the 
help  of  a simple  reciprocity  relation.)  The  system  is  shown 
in  the  accompanying  figure. 
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The  solution  will  be  obtained  as  a sum  of  normal  modes. 
('Ey  a mode  is  meant  a fundamental  solution  of  the  wave  equa- 
tion which  by  itself  satisfies  all  boundary  conditions  except 
those  defining  the  source.)  The  significant  fact  about  this 
form  of  solution  is  that  it  lends  itself  to  practical  numeri- 
cal computation.  In  particular,  no  numerical  integrations 
are  involved.  The  lengthiest  part  of  the  computation,  that 
of  finding  the  phase  velocities  of  the  modes  (i.e.  - finding 
the  eigenvalues),  while  often  somewhat  laborious,  is  still 
practical:  As  an  examination  of  the  "recurrence  method"  of 

Part  I will  show,  the  difficulty  of  finding  the  phase  veloci- 
ties does  not  in  general  increase  very  rapidly  with  the 
number  of  layers  (a  notable  exception  being  the  contrast 
between  the  two  and  three  layer  cases).  Thus  a six  layer 
problem  is  not  "tremendously  more  difficult"  than  a three 
layer  problem. 


After  the  phase  velocities 
are  obtained,  a method  is  used 
tations  which  differs  basically 
This  approach  greatly  simplifie 
it  is  expected  that  it  will  be 
tion  problems:  For  example,  it 
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of  the  propagating  modes 
for  finding  the  mode  exci- 
from  the  usual  techniques, 
s the  present  problem,  and 
useful  in  many  other  excita- 
shouI>d  lend  itself  readily 
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of  elastic  media.  It  has  thus  far  been  developed  only  to 
solve  for  the  propagating  modes  (far  field  case).  In 
order  to  demonstrate  the  method  most  cle ar ly.  it  will  first 
be  used  to  obtain  the  excitations  for  the  simpler  case  of 
a layered  rectangular  wave  guide  (Part  II H these  results 
then  being  extended  to  point  source  excitation  in  plane 
parallel  layers  by  an  appropriate  image  technique  (Part  III)* 
(If  simplicity  of  illustration  were  not  an  important  consi- 
deration, the  technique  of  Part  II,  which  is  applied  to  a 
rectangular  wave  guide,  could  instead  be  applied  directly 
to  the  parallel  layer  case,  as  described  in  paragraph  9c 
of  Part  II.) 


The  author  wishes  to  thank  Dr.  Franklyn  Levin  for  his 
careful  reading  of  the  manuscript  and  for  his  valuable  sugges- 
tions. 
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Part  I The  Phase  Velocities  of  the  Modes 


1-  Referring  to  Fig.  I*  the  velocity  potential  of  a mode  (solution 
of  yrt  , where  j denotes  the  layer)  has  the  form 

' f '-j-  dt* 

in  each  layer  as  shown.  It  is  seen,  in  order  to  satisfy  the 
free  surface  boundary  condition,  that  £(  can  contain  no  cosine 
term. 
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The  boundary  conditions  are: 
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(Continuity  of  pressure) 

(Continuity  of  normal 
velocity) 
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Where  C equals  the  phase  velocity,  which  (like  k)  is  the  same  i-n 
all  layers.  This  relation  is  a direct  consequence  of  the  separa- 
tion of  variables  in  the  wave  equation  of  the  particular  layer, 
and  independent  of  any  boundary  conditions.  (Thus  the  same  rela- 
tion was  used  for  the  three  layer  case,  as  given  in  Eqns.  (A141) 
of  Reference  (a),  as  for  the  present  n layer  case.)  If  we  apply 
the  two  boundary  conditions  to  surface  Z (lowest  surface; 
we  obtain: 
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Applying  the  boundary  conditions  to  any  of  the  other  sur- 
faces (except  the  free  surface)  gives: 
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3.  Thus,  by  applying  Eqns.  (1)  and  (2)  to  the  lowest  boundary, 
and  (3)  and  (4)  to  successively  higher  boundaries  till  Z,  , 
we  obtain  the  A's  and  B's  in  all  layers  as  functions  of 

k and  A*.  But  another  requirement  that  k correspond  to  a 
mode  is  that  J3,  = O . ^Since  in  the  top  layer 

contains  no  cosine  term.)  The  k values  for  which  B,-o  are 
in  general  best  found  by  graphing  B,  against  k (the  A^ 
value  has  no  effect  upon  the  roots).  (Graphing  proves  to 
be  the  practical  procedure  even  for  the  three  layer  case.) 

It  is  seen  that  only  real  values  of  k which  satisfy  the  condi- 
tion correspond  to  propagating  modes.  (Note  also  that  A*, 
must  be  negative  and.  imaginary  to  represent  a propagating  mode.) 

4.  Summary  of  Part  I 

a)  Results  Obtained 

The  k-s  01  the  propagating  inodes  and  the  linear  rela- 
tions between  the  A's  and  B?s  of  each  mode  have  been  found. 
The  remaining  basic  problem  (in  connection  with  solving  for 
the  far  field)  is  to  find  the  point  source  excitations  of 
these  modes. 

b)  Essential  Principle  Employed 

Starting  at  the  lowest  surface  and  applying  the 
boundary  conditions  in  turn  to  successively  higher  sur- 
faces, the  fundamental  solution  in  each  layer  is  obtained 
in  terms  of  the  solution  in  the  lowest  medium  (i.e.  - in 
terms  of  any  k and  any  Am.  without  any  of  the  other  A's 
and  B's).  But  the  only  way  the  top  layer  solution  can  be 
thus  related  to  the  bottom  solution  and  at  the  same  time 
satisfy  the  free  surface  boundary  condition  is  by  permitting 
only  certain  discrete  k values.  (That  is  the  k values 
which  satisfy  all  these  boundary  conditions,  and  therefore 
correspond  to  modes,  are  simply  the  values  which  sabisfy 


the  single  aquation  ' c • While  only  a single 

squat  ici!~  n-:-.  sd  be  satisfied,  and  therefore  solution  bv 
graphing  is  generally  feasible,  the  equation  of  course 
coes  become  more  complicated  as  the  number  of  layers 
increases  . ) 

5o  Additional  Comments 

a)  A necessary  condition  for  the  existence  of  propagating 
modes  is  that  the  wave  velocities  C>  must  not  be  mona- 
tonice L*y  increasing  with  successively  higher  layers. 

(Or  e_se  there  would  be  no  total  reflection  mechanism 
tc  keep  energy  from  escaping  through  the  bottom.) 

b)  In  general  some  of  the  layers  may  have  imaginary 
corresponding  to  propagating  modes  (except  that  (3^  must 
always  be  negative  and  imaginary) . The  mathematics 

above  is  formally  the  same  whether  or  not  certain  of  the  /3a. 
are  imaginary  (i.e.  - in  a layer  where  /3  is  imaginary 
the  sine  and  cosine  terms  become  hyperbolic). 

c)  It  is  interesting  to  compare  the  A y and  B^.  coefficients 
above  with  the  A,  B,  C,  D quantities  in  Eqns.  (A-), 

(A9),  (A30)  of  Reference  (a).  Though  there  is  a.  super- 
ficial resemblence,  it  will  be  seen  that  these  two  sets 
of  quantities  are  basically  different. 


1 


Part  II 


'ti. , 


- 4 - 


Far  Field  Due  to  a Point  Source  in  Layered  Rectangular 
Guide  ' ' 


It  will  be  useful  first  to  consider  the  problem  of  point 
source  mode  excitation  in  the  rectangular  wave  guide  case. 
Not  only  will  the  result  prove  convenient  for  solving  the 
corresponding  parallel  layer  problem,  but  this  simpler  case 
will  also  serve  to  more  clearly  illustrate  a fundamentally 
different  technique  for  finding  Green's  functions. 


Consider 


the  uniforr 
section  in  Fig.  II: 


rectangular  wave  guide  shown  in  cross 


Let  the  width, 


w,  of  the  guide  be  too  small  to 
permit  any  propagating  modes  other 
than  those  which  are  independent 
of  to  exist  (for  the  fixed  fre- 
quency to  which  we  will  consider). 
Assume  further  a large  number, s, 
of  point  sources  all  at  a depth, 
d,  to  exist  as  shown  in  Fig.  III. 
Let  all  be  oscillating  at  the 
same  frequency  6J  , all  be  of 
volume  displacement  SV  (source 
strength  =£V  ) , and  of  the  same 
spatial  phase  relationships  as  an 
m'th  propagating  mode  wave  travel- 
ing in  the  +%  direction.  Thus: 
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Then  for  s sufficiently  large,  we  conclude  from  simple  phase 
considerations.,  that  the  far  field  radiation  produced  by 
these  sources  consists  of  the  m'th  mode  traveling  in  the  -fx 
direction  with  all  other  modes  negligibly  small  in  comparison; 
the  m'th  mode  wave  traveling  in  the-*  direction  is  also 
seen  to  be  negligibly  small. 

(Note : (a)  Strictly  speaking,  to  ensure  that  all  other 

propagating  modes  be  small,  we  may  stipulate  that  the 
sources  be  spaced  uniformly  at  a distance  apart  less 
than  the  shortest  possible  propagating  guide  wave  length 
for  frequency  co  . 

(b)  Actually  the  derivation  to  follow  could  have 
been  accomplished  even  if  only  one  source  existed,  so  that 
the  other  modes  would  not  be  negligibly  small;  but  the 
present  method  is  simpler.) 

3.  Calling  the  velocity  potential  of  this  radiated  m'th  mode  , 
we  assume  in  addition  that  another  wave  fj  . also  of  the  m'th 
mode  and  traveling  in  the  +%>  direction,  has  been  introduced  by 
sources  at  '*--<*>.  ( will  be  used  for  the  layered  rectangu- 

lar wave  guide  to  distinguish  from  the  used  for  the  hori- 
zontally unbounded  parallel  layer  case.)  Further  assume  that 


T 


- 5 - 


is  a much  stronger  signal  than  : Thus  if  is 

"the  ■x  part"  of  the  (Jj'jt  solution  (that  is  = 

Ve®  aL(6x-uttO$  } ? and  Vj$)  is  "the  ^ part"  of  3>x  ? « 

have  \T; ' |»)>7ii  • Also  let  the  sources  producing  the 
wave  begconstrained  to  oscillate  as  previously  described; 
these  oscillations  being  unaffected  by  the  strong 
signal  passing  by»  (The  word  "sources"  will  always  refer 
to  those  producing  the  wave* ) Assume  finally  that  the 
Cfcijr  wave  is  180  out  of  phase  with  the  <£>zwave.  Then  in 
regions  of  large  +•>  the  field  is  given  by  ()r? 


Therefore  since  the  field  approaching  the  source  is 
given  by  amplitude  7)a  , and  the  (far)  field  traveling  away 
is  given  by  the  slightly  smaller  amplitude  \ it 

follows  from  conservation  of  energy  that  the  sources  are 
absorbing  power  (rather  than  losing  it).  Therefore  we  may 
write : 

jftitr  pcr«rg  received  $/  sources]  = [tcmbr  of  % Fiat>J  - jpowtp  of  ( l^-lpil) 


Where  . "by  "2or  field"  is  simply  meant  the  " ^ field" ; the 
" 7?  field"  terminology  being  introduced  for  conciseness, 
since  the  immediately  following  steps  involve  the  7?}s 
rather  than  the  g>  }s  . We  next  introduce  the  evident  rela- 
tion: 

[power  of  ()Va\~\Vi\  fiW)]  _ [j  -IVxifQj] 

[pCWER  CF  FlFljf] 

Wh ere  the  values  of  and  have,  merely  for  concrete- 

ness, been  taken  at  some  arbitrary  fixed  depth  y - a-  . 
Putting  this  relation  into  Eqn.  Ct)  gives:  * 
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Since  7^. (a)  ^(a)  , this  may  be  written: 
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4.  Now  since  the  self  producfjd  pressure  at  a mathematical  point 

source  is  finite  and  continuous  (unlike  the  particle  velocity), 
and  since  !%)>)>  | ;|?r j , it  follows  that  the  pressure  produced 
at  ?ny  of  the  sources  due  to  the  other  sources,  or  the  self 
produced  pressure is  small  compared  to  that  of  the  Vs.  field. 
Further  the  ‘fj ff  and  7 waves , being  180 w out  of  phase,  produce 
the  weakest  possible  resultant:  field;  so  that  from  Eqn.  (5) 
we  conclude  that  the  phasing  of  is  such  that  the  sources 
are  absorbing  the  maximum  power  possible:  Therefore  since  the 
pressure  at  the  sources  is  essentially  that  of  the  Vn_  field, 
and  this  field  is  phased  to  do  a maximum  of  work  on  the 
sources  (the  net  local  pressure  field  works  against  the  "pul- 
sating volume  displacement”  of  the  source) , we  may  write: 
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Where  Ps  is  the  maximum  magnitude  of  pressure  of  the  field 
occurring  at  the  sources,  and  fV  is  the  volume  displacement  of 
each  source.  Hence  Eqn.  (^A)  may  now  be  written: 
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EQN.  C5B) 


Where  ~fi  and  n/^  are  the  instantaneous  real  values  of  % field 
pressure  and  % component  of  particle  velocity  at  any  point. 

5.  We  have  t6t)  . Therefore  since 

Ar„- ~ and  ^a=/^-  ^ (where  /^-density  in  ^ th  rectangular 

wave  guide  layer),  we  haves 
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Where  fti  equals  in  the  j’th  layer;  the  j subscript  merely 
emphasizing  that  ft  is  given  by  different  expressions  in  the 
various  layers. 
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Also,  we  have: 

p = p ( if =fa\%<-4 

1 ' 3t  'HA*  ^ 0r  w* 


Substituting  from  this  relation  and  from  Eqn.  (6)  into  Eqn. 
(5B)  gives: 
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EQN*.  (5TC) 


If  we  choose  cus  d , then 

)=yj{d)=  <2i,  A^ftdL 

Where  (2n  is  the  coefficient  "designating  the  mode  strength" 

(in  the  tor  layer)  of  the  %v(ave.  It  is  important  to  notice 
that  the  set  of  /3}s  (and  $}s  ) which  correspond  to  modes  in  the 
present  rectangular  wave  guide  case  are  Identical  to  those 
applying  in  the  corresponding  parallel  layer  problem  of  Part  I. 
Further  a set~~of  coefficients  and  apply  in  the  present 
case  which  have  precisely  the  same  linear  relationships  as 
the  corresponding  and  fiv  of  i^art  I,  ^his  similarity  occurs 
as  a result  of  the  fact  that  the  "a.  part  of  the  solution"  is 
the  same  (within  a multiplicative  constant)  in  both~  the  parallel 
layer  and  rectangular  wave  guide  cases. 

For  <t=d,  Eqn.  (5C)  may  therefore  be  written: 


I CLi^ftA  1 = - % MN  (5D) 

i*'  h-‘ 

Since  all  the  sources  reinforce  for  the  m;th  mode  in  question, 
we  may  conclude  from  the  linear  superposition  properties  of 
amp  11  txides . that  the  corresponding  expression  for  the  CL,  coeffi- 
cient for  a single  point  source  is: 


BQN  (5’E) 
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Where  AV  is  the  volume  displacement  of  the  source, 
represents  not  the  m'th  mode  produced  by  the  source,  but  a 
Mathematically  fictional*1  m’th  mode  of  arbitrary  strength. 
V*i* fcO  is  the  corresponding  ''fictional  HTCdG'1  Value  at 
2 = d =•  (depth  of  source}.  'Ovi**'  constitutes  a dummy 
v^iable  which  drops  out  (essentially  cancels  the  ??  £7)*  ) 

when  the  evaluation  of  !&>**}  is  carried  through.  ’ 


Thus  Eqn.  (5E)  determines  the  excitation  of  any  propa- 


and  the 

linear  relations 
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the 

and 

j ^ , 

for  the  case 

of  a point  source 

in  the 

top  layer 

of 

an  n 

laver  rectangular  guide.  (Remembering'  that  the  **  are 

identical  with  those  of  the  corresponding  parallel  layer 
problem  of  Part  I,  and  that  the  set  of  (i ^ and  $L,*.  coeffi- 
cients have  the  same  linear  relations  as  the  A and  B 


_ _r»  r>.  ..  i.  -r-  \ 

ui  ret r t ± 


A 


While  the  linear  relations  between  the  Ci^. and  (Bj** coeffi- 
cients are  obtainable  by  the  method  of  Part  I,  and  the 
magnitude  of  £l,^ls  determined  by  Eqn.  (5E),  the  problem  of 
the  relative  phasing  (temporal)  between  the  different 
modes  at  any  value  of  ^ has  yet  to  be  solved:  To  do  this, 

we  observe  that  the  fictitious  0% ^ fields,  corresponding 
to  the  various  propagating  modes  of  a fixed  frequency  (o  , 
all  bear  the  same  phase  relationship  to  the  source  (since 
all  do  maximum  work  on  it).  That  is  all  the  modes 

are  in  phase  with  each  other  at  the  source.  But  since  all 
the  radiated  propagating  modes  are  in  opposite  phase  to 
their  corresponding  fictitious  fields,  these  radiated 

modes  must  also  all  be  in  phase  at  the  source.  Therefore 
w e have  for  any  propagating  mode : 

-^n  w.*  fcff  y awfla tg  /^*/n  ^ 7-  — w ( -«*>. 

Where  is  a constant  which  describes  the  phasing  of  the 
source.  Thus  we  may  write  the  solution: 


(7) 

The  phase  angle  of  ^may  be  completely  accounted  for  by 
the  sin(&w>-tvtvPcj  factor;  so  that  the  coefficient  (%*». may 
be  taken  as  real.  And  from  the  relations  given  in  Part  I, 
being  real  causes  all  the  other  J^»and  fi^to  be  real. 
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8 . Summary  of  Part  II 


Q J 


for  pi? for 


point  source  excitation  of  the  n layered  rectangular 
wave  guide  shown  in  Fig.  II.  More  specifically,  we 
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have^obtained  the  CiJm.  coefficient j so  that  all  the  Cifm. 
and  may  be  obtained  by  the  method  of  Part  I.  Also,  the 
relative  phasing  of  the  different  propagating  modes  has 
been  obtained. 


b)  Essential  Principles  Employed  in  the  Analysis:  The  net 

energy  radiated  or  absorbed  by  a source  is  seen  to  be  the 
work  done  by  or  on  the  source  as  it  ’’pulsates  against  the 
surrounding  pressure  field."  Now,  while  the  work  of  the 
source  against  its  self-produced  pressure  field  cannot  be 
calculated  (the  field  produced  by  the  source  being  pre- 
cisely the  unknown  we  wish  to  find),  the  work  done  when 
a very  strong  and  known  externally  Introduced  signal  travels 
past  the  source  can  to  good  accuracy  be  calculated:  The 

unknown  self- produced  field  of  the  source  contributing 
negligible  fraction  of  the  work  in  the  latter  case. 


4*  Vi  4 r* 
unx  wJ 


fho 


at 

If 
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the  point  source  is  finite,  unlike  the  particle  velocity.) 
the  pulsating  source  is  phased  so  as  to  absorb  energy  from 


by  the  source  must  combine  with 


4- j~\ 


wave  then,  by  conservation  of  energy,  the  field 

externally 
wave  traveling 
incident  signal. 

source  by 


net 


the  oncoming 
radiated 

produced  wave  so  that  the  power  of  the 
away  is  less  than  that  of  the  original 
Thus  by  knowing  the  work  done  on  the  pulsating 
a suitable  strong  field  passing  over  it.  the  field 
by  the  source  is  determined.  A more  precise 


the  principles  involved  was  given  in 


radiated 

description  of 
paragraph  3. 


9 • Additional  Comments 


a)  While  the  case  treated  was  for  the  source  in  the  top  layer 
of  the  wave  guide,  it  is  seen  that  the  same  technique  could 
be  applied  for  the  source  in  any  layer. 

b)  It  is  seen  that  only  the  propagating  modes  are  obtained  by 
the  method  par  se.  That  is,  only  the  far  field  solution 
is  obtained.  No  attempt  has  as  yet  been  made  to  obtain 
either  the  non-propagating  mode  excitations  or  what  corres- 
ponds to  the  branch  line  integral  appearing  in  Eqn.  (A46) 
of  reference  (a) . 

c)  The  technique  employed  here  could  have  been  applied  directly 
to  the  parallel  layer  case  by  considering  the  effect  of  an 
externally  oroduced  cylindrical  mode  in  the  form  of  an 
incoming  wave  as  it  ’’collapses  on  a source"  which  is  con- 
strained to  oscillate  in  a prescribed  manner.  However, 
since  clarity  of  illustration  is  important,  the  method 

was  applied  to  the  considerably  simpler  rectangular  wave 
guide  case.  (For  example,  the  rectangular  wave  guide  case 
was  simpler  in  that  it  was  easier  to  show  that  each  mode 
can  be  "considered  independently"  from  the  others.) 
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Part  III  Far  Field  In  n Parallel  Liquid  Layers  Due  to  Point 
Source  Excitation. 


1.  In  Part  II  we  found  the  point  source  excitations  of  the 
propagating  modes  for  the  case  of  layered  rectangular 
guide.  This  result  will  now  be  extended  to  the  case  of 
point  source  excitation  in  n parallel  liquid  layers.  To 
relate  these  two  cases  we  first  consider  the  following 
image  equality: 
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This  image  equality  states  that  the  field  seen  in  the 
rectangular  wave  guide  due  to  the  line  segment  source  is 
Identical  to  the  field  seen  in  a parallel  layer  region 
due"  to  an  infinite  line  source,  as  shewn. 


t 


2.  The  fundamental  m'th  mode  solution  of  the  wave  equation  in 
cylindrical  coordinates  may,  for  the  top  layer,  he  ex- 
pressed in  the  real  form  (which  will  be  somewhat  more  con- 
venient in  what  follows): 

= A tm,  Jo  (*Ln)  ( ^(wt  *♦ 

Here  A.*,  is  taPen  as  real,  the  phase  angle  Sm.  11  absorbing 
any  phase  angle  that  A*,  may  otherwise  have".  (Ultimately, 

- 6„ ^ will  depend  upon  the  phasing  of  a point  source  with 
respect  to  the  time  t.)  A consequence  of  A/<*being  real 
is  that  all  the  Ai~andB^~  of  the  mode  will  be  real:  That 

is,  all  the  ^»ofrthe  mode  are  real.  The  subscript  on 
Indicates  the  possibility  that  -6~  depends  upon  the  mode. 
(Biit  it  will  turn  out  to  be  independent  of  mode.)  It  is 
seen  that  the  above  expression  for  represents  a standing 
wave.  Using  the  asymptotic  form  of  Jl(A/y  for  w large,  the 
corresponding  far  field  solution  becomes: 


6 = A, /»*■ /lyl  1 c «(ot+sj 

• I AW  tfrx  r i y si  !/  ! * 


The  outgoing  wave  part  of  the  far  field  standing  wave  is: 
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= A EQNU) 
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It  is  this  outgoing  wave  part  that  corresponds  to  radia- 
tion from  a point  source,  the  incoming  wave  part  being 
discarded  (as  is  usually  done  with  the  "physically  unrea- 
sonable" advanced  potential).  (Note  that  A((W,and  are 
the  symbols  used  in  the  cylindrical  coordinate  case  as 
against  the  corresponding  (£,*,  and  j£l<m  of  the  rectangular  case 
of  Part  II.)  Referring  now  to  the  equivalent  parallel 
layer  image  of  Fig.  IV,  the  net  field  produced  (by  all  the 
modes)  at  a large  distance  b from  the  line  source,  and  in 
the  top  layer,  is; 

prom  vEionry  menTiALls.S'=  Y"  ( l <U 

L J 4_u«7-, f^nr  , v 


ALL  MCUS 


( dAJi 


J 


Y/here  ( is  a constant  for  a given  mode  and  is  analagous 

to  the  Aim  coefficient  of  the  "cylindrical"  case.  More 
precisely  the  effective  Aim.  coefficient  of  a short  element 
of  the  line  source  ( Sfr  acts  essentially  like  a point 
source)  would  have  the  value  ^ ^ . 


(9) 


From  the  image  equality  of  Fig.  IV  the  quantity J of  Eqn 
can  be  seen  to  represent  the  total  velocity  potential  inside 
of  a rectangular  wave  guide  (due  to  a line  segment  source): 
This  physical  fact  can  lead  to  evaluating  the  integral  in 
Eqn.  (9)»  (Further,  as  is  often  the  case  in  optics,  this 
integral  has  the  simplifying  property  of  being  evaluable 
over  a stationary  phase  region  given  by  Jpt-)  4c  Jt-  , More 
precisely,  calling  the  "contributing  range  of  v"  — ^ ^ , we 
have  as  jjr->oo  that  ^->00  , but  so  that  fttac.  -+0 . ) 

However,  the  integral  has  essentially  been  evaluated  by  a 
relation  given  on  pg.  416  of  Y/atson’s  "Bessel  Functions", 

2nd  Edition.  The  equation  is: 


f lM)i 
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If  we  put  gl-0 

and  real  for 
of  this  equation  becomes: 


, the  integrand  becomes  imaginary  for 

so  that  for  the  a.*o  case  the  real  part 
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vor  b large  "his  may  be  writ  tent 
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If  this  last  equation  is  expressed  in  complex  form  and 
multiplied  by  Q -U<->t+  ? the  real  part  of  the  resulting 

equation  is: 


V lrt*(  & T I »y+ALx  — 2-  ""U)t  ■€— ) 

- - \ • » » ✓«/  f Y>  Af-  T 

V 


EQhi  (iO) 


Identifying  this  integral  with  that  appearing  in  Eqn.  (9), 
gives  forj'  : 

•J  = V A *£&&..  (j*ULlr-at-<b 
H'  iL 


Since  « is  also  the  field  seen  in  the  rectangular  wave  guide 
(by  the  imago  equality  of  Fig.  IV) j and  since  the  summand  is 
of  the  form  fccuauif] . c^. OL,l-~vt ~ ^0  ? the  summand  is  recognized 
as  constituting  the  m!th  rectangular  wave  guide  mode*  (i.e., 
The  summand  is  an  undamped  sinusoidal  function  of  , as 

is  the  case  for  propagating  rectangular  wave  guide  modes.) 
Thus  we  may  writes 


jP^ecfA^GOiAR  WWt  QWE  VZi-OCf T/  POTEMTIU 
Or  /m‘tb  Mcpe  AT  % -Jr  4*5  sN  Tc?  l/.JOR  J 


* ^ = (IfKx)  ^L> 


Next  we  observe  that  may  be  regarded  as  the  rectangular 
wave  guide  m;th  mode  solution  for  a point  source  of  a volume 
displacement  £v  equal  to  that  of  the  segment  source  of  Fig.  IV 
(of  length  w).  (The  far  field  equivalence  for  a rectangular 
wave  guide  between  the  line  segment  source  and  a point  source 
of  equal  displacement  can  be  understood  from  the  fact  that  the 
same  work  is  dons  on  either  type  source  by  the  "hypothetical 
field"  of  Part  II,  provided  only  that  w is  small 
enough  so  that  the  propagatin,  m ,de  amplitudes  are  independent 


- 13  - 


of  y.  (The  value  of  w in  no  way  affects  the  equivalent 
parallel  layer  image;)}  Therefore  we  have  that  ■&,  may 
also  be  expressed  by  Eqn.  (7)* 


Where  (2/*.is  given  by  Eqn,  (5E)  when  in  that  equation  is 
put  equal  to  Xtjp  > "here  (&)  ls  the  strength  per  unit 

length  of  the  line  source;  and  where  H is  some  phase  angle 
indicating  the  fact  that  the  phasing  of  the  line  source  has 
not  been  defined.  Comparing  the  coefficients  of  these  last 
two  expressions  for  ^ gives: 


EQN  (") 


Eqn.  (11)  relates  the  unknown  (zd"/1  of  the  cylindrical  case 
to  the  known  CL**,  of  the  rectangular  wave  guide  case.  It  " 
leads  to  a connection  of  the  point  source  mode  excitation 
in  the  cylindrical  case  to  the  previously  determined  excita- 
tion in  the  rectangular  guide  case. 


5»  Eqn.  (11)  may  be  written: 


dAi 
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'Where  the  volume  displacement  per  unit  length  of  the 

Dividing 

, _ \ elves : 


line  source  of  Fig.  IV  (as  already  defined), 
through  by (4^  ) gives : 


$ A irm 

w 


Where  is  the  effective  Am*  coefficient  (cylindrical 

case)  due  to  radiation  from  the  differential  line  source 
element  SV,  And  where  SV  is  the  total  volume  displacement 
of  the  element  » Identifying  the  differential  element 
St*  with  an  effective  point  source,  and  calling  STW^Am* 
(to*  correspond  to  the  nomenclature  for  a point  source  in 
the  cylindrical  case),  the  last  equation  may  be  written: 
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To  obtain  \a  !* I 
in  Eqn.  (5E) : 


I (CL  is  real)  we  re  place  A\J  by 
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Putting  this  value  for  JCLi^l  into  Eqn*  (11A)  gives: 


Aim 


%~dt 

lh 


EQN  (M  b) 


Eqn-  (11B)  gives  the  magnitude  of  /A.—  for  a point  source 
of  volume  displacement  SV  and  at  a depth  d in  the  top 
layer  in  the  parallel  layer  (l.e.  - cylindrical)  case- 
The  other  quantities  appearing  in  Eqn.  (11B)  are  identical 
with  those  defined  in  connection  with  Eqn.  (5B). 


6.  The  temporal  phasing  of  the  modes  due  to  point  source 

excitation  in  the  parallel  layer  case  is  yet  to  be  deter- 
mined: 

Comparing  the  integrand  to  the  value  of  the  integral  in 
Eqn.  (10),  we  see  that  the  signal  from  a line  source  a 
distance  b away  (corresponds  to  integral)  lags  the  signal 
from  a point  source  the  same  distance  away  (corresponds 
to  the  integrand  when  y=o)  by  Vf  Tn  any  given  mode. 

And  by  the  image  equality  of  Fig.  IV  the  mode  phasing  from 
the  infinite  line  source  is  the  same  as  for  the  rectangu- 
lar wave  guide  case.  Thus  the  far  field  m‘th  mode  signal 
due  to  a point  source  in  the  parallel  layer  case  is  of 
phase  angle  ^ less  than  the  corresponding  (same  distance 
from  similarly  phased  source)  signal  in  the  rectangular 
guide  case. 


7.  Let  us  assume  that  the  point  source  is  phased  in  time  so 
as  to  have  its  instantaneous  volume  displacement  varying 
as  ccn  oj t . First  we  consider  the  rectangular  wave  guide 
case  for  this  source:  Referring  to  Part  II,,  since  the 

field  does  a maximum  of  work  on  the  source,  its  pressure 
can  be  seen  to  lag  the  instantaneous  source  displacement 
by  90° » Therefore  the  radiated  mode,  being  out  of  phase 
with  the  fie Id , must  have  its  pressure  leading  the 

source  displacement  by  9CU • And  since  , It 

readily  follows  that  the  instantaneous  soi  rce  displacement 
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is  in  phase  with  the  velocity  potential  (at  the  source) 
of  any  of  the  radiating  modes.  Hence  the  temporal  phase 
of  the  velocity  potential  of  a mode  due  to  the  source  may 
be  expressed  at  a point  (%%)  in  the  top  layer  by* 

Hie 

The  ± factor  indicates  that  an  additional  180°  phase 
shift  between  observer  and  source  must  be  introduced 
when  their  respective  depths  are  such  that  and 

are  opposite  in  sign.  (Since  the  temporal  phas- 
ing of  a given  mode  alternates  with  depth  as  the  sign  of 
Further,  as  is  required,  the  argument  of  the 
cosine  is  such  that  if  the  observer  is  at  the  source 
(.%-o ) 2-*-  d)  the  expression  becomes  identical  to  that  for 
the  phase  of  the  instantaneous  source  displacement: 

<p  (Qjdjt)  + Zo± (-wtr)  ==  r C 

/7*\ 


8. 


From  paragraph  6 it  follows  we  may  obtain  the  correspond- 
ing parallel  layer  phasing  (far  field)  simply  by  sub- 
tracting from  the  phase  of  ^ ( and  replacing  the  'X- 
coordinate  by  A )s 


/'JfctFef?  y-tctufi 


Thus  we  know  the  phase  of  a far  field  radiated  mode  in 
the  parallel  layer  case,  and  from  Eqn.  (11B)  we  know  the 
magnitude  of  the  corresponding  Ai/m.  coefficient.  These 
results  may  be  expressed  more  conveniently  if  we  generalize 
Eqn.  (11B)  by  defining  A^  as  a (real)  quantity  whose  sign 
depends  upon  the  depth  of  the  source: 


EQrt(n 


Where  this  value  of  Aim  is  to  be  inserted  into  the  follow- 
ing solution  for  the  far  field  radiated  mode  in  the  top 
layers 


r^\rr\ 

[V<<> 

A->oo  ’ V ' #Vw-  / 
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The  form  of  this  last  equation  can  be  seen  to  follow 
directly  from  the  corresponding  expression  in  paragraph  2 
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In  general  the  solution  in  any  of  the  layers  is  given  by: 


£QW(I3) 


Where  the  are  defined  in  Pig.  I (in  Part  I).  (The  m 
subscripts , representing  the  mode,  do  not  appear  in  Fig.  I.) 


(p.  = £ 
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9 « Summary  of  Part  III 

a)  Results  Obtained:  Solved  the  far  field  problem  for 

point  source  excitation  for  the  case  of  n parallel 
liquid  layers.  More  specifically,  -we  have  obtained 
the  Aim,  coefficient  (defined  in  Fig.  I)  for  the  case 
of  a point  source;  so  that  all  the  and  B*m-  may  be 
obtained  by  the  method  of  Part  I.  Also  the  phasing  of 
the  propagating  modes  with  respect  to  the  phasing  of 
the  source  have  been  obtained. 

10  b)  Essential  Principles  Employed  s It  will  be  simpler  if 

we  describe  a slight  variation  of  the  method  employed. 
(Though  the  method  employed  was  more  convenient  for 
the  actual  derivation.)  Understanding  the  variation 
wiil  basically  amount  to  understanding  the  original 
method.  If  in  Fig.  IV  the  line  source  segment  in  the 
rectangular  wave  guide  is  replaced  by  a point  source 
midway  between  the  guide  walls,  then  the  equivalent 
Image  is  seen  to  consist  of  an  infinite  "linear  string" 
of  point  sources  (uniformly  spaced  at  a distance  apart 
equal  to  the  width  of  the  rectangular  guide).  There- 
fore we  may  connect  the  "cylindrical  solution"  Tun-- 

k now n excitation)  to  the  rectangular  guide  solution 

fe x 1“ r- i - n known  from  Part  II)  as  fellows i 


1 10*1/  i lON 


rca  PC'iN  r SCV3LE  IN  tieCTVKUJiH  QOitjj  — ^^SCUTiQM  FCR  EACH  PONT  jN  PAJtMiet 


Aa  PcfMr 
sovset 


S n,  je  ihe  left  hand  side  Is  known  from  Part  II „ it 
nurs  Out:  she  "cylindrical  mode"  excitations,  which 
a 0£i§i .*■  in  che  right  hand  side  can  be  determined . 


To  show  one  procedure  in  somewhat  more  detail,  v/e 
Write  (he  last  expression  for  the  top  layer  as  folcowsb 
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Upon  evaluating  the  inner  summation  on  the  right  hand 
side  it  is  found  that  the  (remaining)  summand  may  be 
equated  to  the  summand  on  the  left  hand  side,  so  that 
the  Aim  (the  "excitations")  may  be  obtained. 

10 o Additional  Comments 

a)  The  same  comments  appearing  in  paragraphs  6a  and  8b 
of  Part  II  apply  here. 
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RESULTS  AND  CONCLUSIONS 


The  problem  of  the  acoustic  wave  produced  in  N lossless 
plane  parallel  liquid  layers  due  to  point  source  excitation 
has  been  solved  for  the  far  field  in  a form  permitting  of 
practical  computation.  The  solution  is  expressed  as  a sum 
of  normal  modes:  This  constitutes  a generalization  of  the 

results  of  Pekeris  who  obtained  solutions  in  the  normal  mode 
form  for  the  special  cases  of  two  and  three  layers  (Reference 
(a)  ) . 

The  method  employed  here  for  obtaining  the  mode  excita- 
tions was  basically  different  from  the  usual  techniques  of 
finding  Green !s  functions.  This  method  greatly  simplified 
the  present  problem,  and  may  similarly  be  expected  to  facili- 
tate the  procedure  in  various  other  excitation  problems. 
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Appendix  A: 


Computation  Proceduie 


For  practical  computing  purposes,  as  v»ell  as  for  the  sake  of 
conciseness,  the  results  will  now  be  summarized  in  a more 
directly  applicable  form: 

1 . Statement  of  Problem 

To  obtain  the  far  field  velocity  potential  produced 
in  a system  of  n parallel  dissipationless  liquid  layers 
by  a point  source  of  sound  oscillating  at  a frequency  , 
The  system  is  shown  in  the  following  diagram. 


JLh. 

(I,  Ca 


‘Wrcc 


FkSt  SVKFM& 


^ j £>»» 


iem  -iHftMtrr 
liquid  ^ 


A necessary  condition  for  the  existence  of  propagating  modes 
is  that  the  wave  velocities  Oj.  must  not  be  raonatonically 
increasing  with  successively  higher  layers. 


'orm  of  the  Far  Field  Solution 


. = t ( J 

FitLD)  I * 


&*(<■>* 


VJ here  the  ^ subscript  refers  to  the  layer  (as  used  for  the 
/O's  and  in  the  above  diagram)  and  the  m subscript  to 
the  mode. 
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3 . Finding  the  of  the  Propagating  Modes,  and  the  Linear 

Relations  Between  the  Aj*.  and  Bfm*  of  Any  Given  Mode* 

First  obtain  the  as  functions  of  A (eigenvalues 
of  A unknown  as  yet)  and  the  Cj,  using  the  relations 
given  in  the  preceding  paragraph  for  the  : Here  the 

m subscripts  are  omitted  because  the  are  regarded 
temporarily  as  continuous  functions  of  A , rather  than 
as  discrete  values  corresponding  to  modes. 

Then  apply  Eqns . (1)  and  (2)  of  Part  II  to  the 

lowest  boundary  and  Eqns.  (3)  and  (4)  to  each  successive- 

ly higher  boundary  through  . This  gives  the  f\L  and  Bj. 
as  functions  of  A and  A~.  0 

Next  we  must  find  the  lvalues  for  which  B^Aj ArA) = & 
These  A values  correspond  to  the  propagating  modes,  and 
are  designated  by  a subscript  as  The  corresponding 

& , f\j.  anu Dj.  values  are  thus  designated 
respectively.,  If  three  or  more  layers  are  involved,  the 
A values  for  which  B,~0  are  in  general  best  found  by 
graphing  B,  against  & . (The  unknown  value  has  no 
effect  upon  the  roots.)  Only  real  values  of  A which 
satisfy  this  condition  correspond  to  propagating  modes. 

4.  Finding  the  Coefficient  /W  ; And  Consequently  Obtaining 
All  The  and~]fy~  of  The  Mode . 

The  coefficient  A 1/1*.  is  given  by  Eqn.  (11C).  The 
other  A fin.  and  of  the  mode  may  then  be  determined  bv 

virtue  of  their  relations  to  /A)/Wv  through  Eqns.  (1),  (2), 
(3)>  (4) • (Eqns.  (1),  (2),  (3)?  (4)  are  such  as  to  ylelo 
readily  the  /\^m  and  in  terms  of  the  bottom  coeffi- 
cient but  must  be  essentially  "Inverted"  In  order  to 
give  the  A}**  and  Bjm.  in  terms  of  /\,^  instead.  This 
inversion  of  the  equations  is  a very  simple  iterative  pro- 
cess and  will  not  be  done  here.) 


5 • Far  Field  m*th  Mode  Solution  for  Outgoing  Wave  Due  to  a 
Point  Source  in  the  Top  Layer. 

Having  obtained  all  the  fym.  , Afm-  •>  and  correspond- 
ing to  each  ^ we  have  the  solution  for  the  velocity  po- 

in  the  j'th  layer  by  substituting 


we 

tential  of  the  m'th  mode 
into  the  expression  for 
above . 


6.  , 

(Pk'IVOjlfi) 

Kit*  D ) 


given  in  paragraph  2 
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Appendix  B:  Comparison  of  Two  Layer  Case  to  Results  of 

Reference  (sQ.  . - — 


1.  For  the  two  layer  case  (i.e.  - a single  layer  over  a semi- 
infinite region)  Eqns.  (1)  and  (2)  of  Part  I becomes 


£ Qbi  (l  A) 

B,  = 0=  A.e  ^-a 

eqn  (*a) 

From  Eqn.  (2A)  we  have 

or 

A- A?.  ■ ‘1  £ 

eqn  (^b) 

If  the  fis  in  this  last  equation  are  expressed  in  terms  of 
those  values  of  H satisfying  the  relation  correspond  to  the 
propagating  modes.  If  expressed  in  the  appropriate  notation 
this  equation  becomes  identical  with  Eqn.  (A49)  of  reference 
(a) . 


2.  Eqn.  (HC)  of  Part  III  becomes  for  the  two  layer  case  (omit- 
ting the  m subscripts  will  nxst  cause  confusion): 


-o 

CO 


(tZi 


<2.  Aiyxfid. 


Where  di  and  bZaare  coefficients  corresponding  to  a 
“mathematically  fictional11  m'th  mode  of  arbitrary  strength. 
That  is  the~ZZ7  and  #aObey  Eons.  (1A)  and(2A).  The  ex- 
pression readily  becomes: 


CO  £V  /wn/3,c(. 
-v  — aA-  5 > 


j-Hai  e'*m' 

LP‘  a; 
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Since  &i  and  (X2  satisfy  Eqns.  ( Ik ) and  (2A),  we  may  write 
Eqn-  as: 


a, 

a, 


or 


A, 

A, 


! 


hence 


cj  8\J  Aiti (3, d. 


Mri  20' fr 


+ 


_J_ji 

^ rf 


Introducing  the  notation  used  by  Pekeris  (Reference  (a)), 
we  have: 


/I 


Using  Lqn.  (2B)  above,  we  get,  in  this  notation: 


A,  - 


< ’X)  SV  AbnficL 


|-J  — Xjbw  ^-'/ry  J <£/" 

ft)  \ X<* 


which  reduces  to: 


a-  Co-*  lUt 


'X,,  /U/rt  "X.-n.  ^ Xn. Aw? 
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Comparing  to  the  definition  of  F(%*)  given  in  Eqn.  (A73) 
of  Reference  (a),  we  may  write: 


A 


cj  SV 
H 


Therefore  the  solution  as  given  in  paragraph  2 of  Appendix 
A becomes,  for  the  top  layer: 


<P, 


i/m  (otTCi'.VC-  WAvfj 
Far 


. H 


Putting  this  into  complex  form, 
this  becoms: 


and  writing  $ - 


y^ 

~FT  » 


If  the  source  strength  ^ is  put  equal  to  this  solution 

becomes  identical  with  Pekerls1  result  (Eqn.  A7l). 


3.  The  question  of  whether  the  present  solution  is  identi- 

cal with  that  of  Pekeris  depends  upon  whether  nr  not  he 
chose  the  source  strength  as  , as  is  required  by  the 
above  comparison  of  results.  At  first  sight  the  two  results 
would  appear  contradictory  sines  does  not  have  the 
dimensions  of  volume  appropriate  to  a source  displacement* 
However  there  is  actually  no  disagreement,  for  it  will  be 
shown  that  Pekerls  nevertheless  did  in  effect  fi.e.  - im- 
plicitly) choose  the  source  strength  as  ttgj  : Since  there 

was  no  particularly  important  reason  to  do  so  at  the  time, 
he  did  not  attempt  to  keep  his  equations  dimensionally 
consistent;  and  by  introducing  the  source  strength  indirectly 
through  a boundary  condition,  an  implicit  consequence  was 
that  the  source  strength  was  #£■  . To  show  this  we  proceed 

as  follows: 


We  will  make  use  of  equations  (A17)  and  (Al6)  of 
Reference  (a).  We  notice  that  each  of  these  is  dimension- 
ally  inconsistent,  its  right  hand  side  having  the  dimensions 
of  length  while  the  left  hand  side  has  the  dimensions  of  a 
velocity  potential  (- f)  . The  discontinuity  in  the 
particle  velocity  at  the  plane  (which  is  a olaitff 

containing  the  source)  has  a vertical  component  given  by 


xatUJft  ' _ a; 
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Using  the  values  of  Vi  and  i'-A  given  by  equations  (Al/ j and 
(A18)  this  discontinuity  is  C'eadily  shown  to  be  equal  to 
j"  r{ft*c.  Rctk  . This  ®-s  a discontinuous  integral, 

being  infinite  at  and  z®ro  everywhere  else.  However, 

if  this  discontinuous  integra^  is  integrated  over  the 
plane  the  result  is  finite.  ^ this  multiple  integration 
is  carried  out  (for  example  with  the  help  of  equation  (6), 
pgo  386  of  Watson1  s **Bessel  Sanctions11 , 2nd  Edition)  the 
value  is  found  to  be  ^JT  • Knowing  the  integrated 
value  of  the  vertical  velocity  discontinuity  over  the  plane 
'$■  - d-  s we  can,  bv  simple  flc1^  c. ,9.U s.  1r ,d,§.r_%t i-Q ns , show  that  the 

maximum  volume  displacement  ^ 

tty  — 2 ji_ . o However  equati“^8  (A17)  and  (A18)  represent 

the  standing  wave  solution.  ^ only  the  outgoing  wave 
part  of  the  solution  is  kept  i^  can  shown,  as  may  be  ex- 
pected, that  the  source  strength  becomes  halved?  ^ . 

But  this  Ls  the  value  that  vas  required  to  verify  the 
equivalence  between  the  pres0611^  solution  and  that  of  Pekeris. 


Thus,  we  may  conclude  t^a^  the  present  results  are  in 
agreement  with  those  of  Pekeuus. 
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